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Abstract 

We unfold the codimension-two simultaneous occurrence of a border-collision bifur- 
cation and a period-doubling bifurcation for a general piecewise-smooth, continuous 
map. We find that, with sufficient non-degeneracy conditions, a locus of period- 
doubling bifurcations emanates non-tangentially from a locus of border-collision bi- 
furcations. The corresponding period-doubled solution undergoes a border-collision 
bifurcation along a curve emanating from the codimension-two point and tangent to 
the period-doubling locus here. In the case that the map is one-dimensional local 
dynamics are completely classified; in particular, we give conditions that ensure chaos. 

Piecewise-smooth, continuous maps are used to model nonsmooth physical 
situations with discrete-time inputs, and also arise as Poincare maps in a variety 
of piecewise-smooth, systems of ordinary differential equations. Frequently a 
period-doubling bifurcation is the cause for a critical change in dynamics of a map 
producing a physically important period-doubled solution. A special situation 
arises in a piecewise-smooth, continuous map if this bifurcation occurs at a non- 
differentiable point in phase space. We show that soon after its creation, the 
period-doubled solution undergoes a border-collision bifurcation. Beyond this 
bifurcation the period-doubled solution may persist but the solution may also 
immediately become chaotic. 

1 Introduction 

Piecewise-smooth systems are currently being utilized in a wide variety of fields to model 
physical systems involving a discontinuity or sudden change. Examples include electrical 
circuits with a switching component such as a diode or a transistor [HE113], vibro-impacting 
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systems and systems with friction [U [5j El [7J [8] , and optimization in economics [9j [TOj [TT] . In 
this paper we study discrete-time, piecewise-smooth systems that are everywhere continuous. 



is piecewise-smooth continuous if F : M — > R is everywhere continuous but non-differentiable 
on co dimension-one boundaries called switching manifolds. The collision of an invariant set 
with a switching manifold under parameter variation of (pQ) may produce a bifurcation that 
cannot occur in smooth systems. When the invariant set is a fixed point and a bifurcation 
occurs, the bifurcation is called a border- collision bifurcation. 

Border-collision bifurcations may be analogous to familiar smooth bifurcations, such as 
saddle-node and period-doubling bifurcations. Alternatively they may be very complex. For 
instance a chaotic attractor may be generated at a border-collision bifurcation, even if the 
map is only one-dimensional [12] . In two or more dimensions, a complete classification of all 
possible dynamical behavior local to a border-collision bifurcation is yet to be determined, 
see for instance [131 HI H51 ESI Ej . 

At a border-collision bifurcation a fixed point, x*, lies on a switching manifold. Due to 
the lack of differentiability the stability multipliers associated with x* are not well-defined. 
However, the switching manifold locally divides phase space into regions where the map 
has different smooth components; by continuity x* is a fixed point of each component and, 
restricted to each component, associated multipliers of x* are well-defined. Typically x* 
will be a hyperbolic fixed point of each component. A special situation arises when x* is 
non- hyperbolic in at least one component. 

Under the assumption that the border-collision bifurcation occurs at a smooth point on 
a single switching manifold, there are exactly two smooth map components (half-maps), in 
some neighborhood. Consider the situation that x* is hyperbolic in one half-map but non- 
hyperbolic in the other half-map. Then there are three co dimension-two cases to consider; 
namely that on the non-hyperbolic side, x* has an associated multiplier 1, a multiplier 
— 1 or a complex conjugate pair of multipliers, e ±2nmj . These correspond to the simultaneous 
occurrence of a border-collision bifurcation with a saddle-node bifurcation, a period-doubling 
bifurcation and a Neimark-Sacker bifurcation, respectively. Unlike a generic border-collision 
bifurcation, in each codimension-two case linear terms of the half-maps are insufficient to 
describe all local dynamical behavior. 

For the first case, generically a locus of saddle-node bifurcations emanates from the 
codimension-two, border-collision bifurcation which in a two-parameter bifurcation diagram 
is a curve tangent to a locus of border-collision bifurcations [18] (see [19J for the continuous- 
time analogue). The third case may be very complicated and is left for future work. The 
precise nature of local bifurcations may depend on the rationality of to. A recent analysis 
of these situations in a general piecewise-smooth map is given in [2U]. The purpose of the 
present paper is to determine generic dynamical behavior of (pQ) local to the codimension-two 
border-collision bifurcation of the second case. 

The coincidence of border-collision and period-doubling has been observed in real-world 
systems. Recall that a corner collision in a Filippov system provides an example of a border- 
collision bifurcation in a map [21J. In [221 123] the authors describe the coincidence of a 
corner collision with period-doubling in a Filippov model of a DC/DC power converter. 



A map 




(1) 
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Grazing-sliding (which also corresponds to a piecewise-smooth, continuous Poincare map, 
see (21]) and period-doubling have been seen to occur simultaneously in a model of a forced, 
dry- friction oscillator [25] . 

In the next section we introduce a piecewise form for ([!]) suitable for the ensuing analysis 
and compute its fixed points. In §3] we unfold the codimension-two situation for the case 
that the map is one- dimensional. Here all local dynamical behavior is determined. We find 
there are essentially six distinct, generic, unfolding scenarios though each exhibit the same 
basic bifurcation curves. When the map is of a dimension greater than one we prove that 
the same basic bifurcation curves exist but do not classify all possible local dynamics, §HJ 
Finally §0 gives conclusions. 



2 Formulization of the codimension-two point and setup 

We wish to explore local dynamical behavior near an arbitrary border-collision bifurcation of 
(Cp). We restrict ourselves to a neighborhood in both phase space and parameter space of the 
border-collision bifurcation within which we assume the border-collision bifurcation occurs 
on an isolated C l switching manifold and that away from the switching manifold, F is C k . 
Then, via coordinate transformations similar to those given in [SH], we may assume the 
border-collision bifurcation occurs at the origin, x — 0, when a parameter, is zero, and 
that to order I the switching manifold is simply the plane eJx = 0. In this paper we are 
not concerned with effects due to a nonsmooth switching manifold (for studies of various 
piecewise-smooth systems involving a non-differentiable switching manifold, see, for instance, 
[21 | [271 1281 129]). For this reason we assume / is sufficiently large to not affect local dynamics 
and for simplicity assume the switching manifold is exactly the plane eJx = 0. 

For convenience we denote the first component of the iV-dimensional vector, x, by s, i.e. 

s = eJx . 

Let rj be a second map parameter, independent of ji. Then we are interested in local dynamics 
of 

where 

f {i) (x; /i, r]) = fib(fi, 7]) + Ai((i, r))x + 0(|x| 2 ) + o(k) , (3) 

is C k for i = L,R. Here b G M. N and Al and Ar are NxN matrices which, by the assumption 
of continuity of ([2]), are identical in their last — 1 columns. Throughout this paper we use 
the notation O(k) [o(k)} to denote terms that are order k or larger [larger than order k] in 
all variables and parameters of a given expression. We refer to as the left half-map and 
f( R> ) as the right half- map. 

Let us first determine fixed points of (121) . Assume for the moment that for some rj, 1 
is not an eigenvalue of each Ai(0, rj). Then near (x; ji) = (0;0), each half-map, /W, has a 
unique fixed point 

x*« (fj,, ri) = (I- A(0, r ? ))- 1 6(0, V )ji + 0{ji 2 ) . (4) 
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Recall that the adjugate of any matrix X, adj(X), obeys a.dj(X)X = det(X)I [301 El] • Since 
Al and Ar can differ in only the first column, it follows that the first row of adj(J — Al) 
and of adj(J — A R ) are identical. We denote this row by g T : 

g T = e]adj (/ - A L ) = e[adj (/ - A R ) . (5) 

Multiplication of (jlj) by ej on the left yields 



detU - A ) 



M + 0(/i 2 ), (6) 

fi=0 



X 



where s*w denotes the first component of 

The point x*( L > is a fixed point of the piecewise-smooth map (T5J) and said to be admissible 
(that is, a fixed point of (J2])) whenever s*^ < 0, otherwise it is said to be virtual. Similarly 
x*^ is admissible exactly when s*^ > 0. The condition 

Q T bU =0 ± , (7) 

ensures that the parameter ji "unfolds" the border-collision; we will assume this nondegen- 
eracy condition holds. 

Let us briefly review Feigin's classification of border-collision bifurcations by the relative 
admissibility of fixed points and two-cycles (refer to [22] for further details). If / — A L {0, rj) 
and I — Ar(0, rj) are both nonsingular and ([7]) is satisfied, for small ji each x*^ is admissible 
when ji = and for one sign of ji. Let erf [a^] denote the number of real eigenvalues of 
A(0, rj) that are greater than 1 [less than —1]. By ([6]), if a\ + cr^ is even, x*^ and x*^ are 
admissible for different signs of ji. In this case as ji is varied through zero a single fixed point 
effectively persists. If instead cr^ + a^ is odd, x*^ and x*^ are admissible for the same sign 
of /jl; at fj, = the two fixed points collide and annihilate in a nonsmooth fold bifurcation. 
Furthermore if / — A^O, ?7)Ar(0, 77) is nonsingular, —1 is not an eigenvalue of ^(O,^) and 
Ar(0, rj), and + is odd, then a two-cycle of ([2]) (that collapses to the origin as ji — ► 0) 
is admissible for one sign of ji, whereas if a L + cr^ is even, no such two-cycle exists. 

When ji = 0, the origin is a fixed point of (J2J) that lies on the switching manifold, s = 0. 
For small ji 7^ 0, stability of, and dynamics local to, the fixed point, x*^\ are determined 
by the eigenvalues of Ai(0,rj). In nondegenerate situations the eigenvalues of the matrices 
A; (0,7/) completely determine dynamics local to the border-collision bifurcation at ji — 0. A 
special situation arises whenever either ^4^(0, rj) or A R (0, rj), has an eigenvalue that lies on the 
unit circle. Without loss of generality we may assume the former matrix has an eigenvalue 
on the unit circle when 7/ = 0. In this situation a local smooth bifurcation may occur for the 
left half-map, f*^ L \ at ji = rj = 0. The nature of the smooth bifurcation is determined by 
nonlinear terms, consequently such terms are important to the border-collision bifurcation. 

The remainder of this paper is an analysis of the case that —1 is an eigenvalue of Al(0, 0). 



3 Unfolding in one dimension 

We begin by treating the one- dimensional case, that is, x = s £ R. Then the map ([2]) may 
be written as 

_ / f (L) {x n ;ji,rj), x n <0 , , 

n+1 " I f iR) (x n ,ji,rj), x n >0 ■ { * } 
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We assume that the two smooth components of the map are at least C 3 (i.e. k > 3) and 
write them as 

f {L \x; fi, r]) = fib(fi : rf) + a L (n, T])x + p(/j,, rj)x 2 + q(fi, rj)x 3 + o(x 3 ) , 
f {R) (x; /I, rj) = nb(fx, rj) + a R (fi, rj)x + 0{x 2 ) . 

The map (JSj) has a border-collision bifurcation at 11 = 0. We assume that the codimension- 
two situation of interest occurs at // = r\ = 0, that is Ol(0, 0) = —1. In order for parameters 
\x and ?y to unfold the bifurcation therefore we require 6(0, 0)^0 (by©) and ^(0,0)^0. 
For simplicity we assume appropriate scalings upon ii and rj have been performed such that 
these values are both 1 as stated in the following theorem. 

Theorem 1. 

Consider ^Bj) with (TJ|) and suppose 

i) Ol(0,0) = —1 (singularity condition), 

ii) 6(0,0) = 1 (border- collision non- degeneracy condition), 

Hi) ^-(0,0) = 1 (transversality condition). 
Let 

a { R) = a fl (0,0), (10) 
c = p 2 (0,0) + g(0,0) . (11) 

Then there exists S > such that the left half-map restricted to the neighborhood 

M = |(x;/i,r/) |x|, |//|, \rj\ < Sj (12) 

has a unique fixed point given by a C k function, x*( L \fi,r]) = |/i + 0(2), and there exist 
unique, C k ~ l functions hi,h 2 : R — > R £/ia£ satisfy 



/i + 0(// 2 ), (13) 

//=7?=0 

fcaC/i) = - j/i 2 + o(/z 2 ) , (14) 



such that 



1) when ji < ; hi(fj)) has an associated multiplier of — 1, and if c ^ £/ie curve 
77 = /ii(/i) corresponds to a locus of admissible, period- doubling bifurcations of x*( L \ 

2) the origin belongs to a period-two orbit of f^ on rj = h 2 (fi) that is admissible when 

fJL < 0, 

3) if Oq 7^ ±1 and cq ^ ; then fixed points and two-cycles of (EJ) exist in the sectors 
shown in Fig. [IJ 
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Figure 1: Sketches of bifurcation sets of (|HJ) for different values of a and Cq, see Theorem 
[Q Blue curves correspond to loci of period-doubling bifurcations; black curves correspond 
to border-collisions and are shown dashed when no topological change occurs at the border- 
collision. Insets show phase portraits with circles used to denote fixed points and triangles 
used to denote points on two-cycles. These are blue when the solutions are stable and red 
otherwise. The small vertical lines denote the switching manifold. Chaotic dynamics occurs 
only for parameter values in the shaded regions; indeed, periodic solutions of a period higher 
than two only occur in these regions. The bifurcation diagram shown in Fig. [2] is taken along 
the dash-dot line segment of panel F. 
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4) if % < 1 or fi > or r\ > h 2 (fi), (EP has no n-cycles for any n > 3 in N , 

5) if Oq > 1 and fx < and rj < ft^O"), then (Gj) exhibits chaos in M . 



A proof of Theorem [T] is given in Appendix [A] To facilitate the extension of Theorem 
CD to higher dimensions in §H the theorem is stated in a manner that yields a partial result 
when the non-degeneracy condition on the nonlinear terms, Cq ^ 0, is not satisfied. The sign 
of Co determines the criticality of the period-doubling bifurcations. 

Fig. [1] shows the six different bifurcation scenarios depending on the value of and the 
sign of Co that are predicted by the theorem. In each case the nature of the border-collision 
bifurcation at fi = for small rj ^ may be determined by refering to the well-understood, 
one-dimensional, piecewise-linear map 

J ~\~ &LXm X n ^0 ,^ _n 

Xn+l ~ \ V + a R x n , x n >0 ' 

see for instance [12]. For example, along the positive r^-axis in Fig. [T]-A, we have ol ~ 
— 1 + rj > — 1 and or < —1, thus by Feigin's classification a single fixed point persists and a 
two-cycle is created that coexists with the left fixed point and is stable. 

In panels E and F, along the negative 77- axis we have ai ~ — 1 H- 77 < —1 and or > 1. 
At points just to the left of the negative 7/-axis, both fixed points and a two-cycle are 




Figure 2: A bifurcation diagram for ([8]) when b — 1, ul — f] — 1, Or = §, V — —1, 1 = \ an d 
there are no additional higher order terms. The value of r\ is fixed at —0.25. The parameters 
considered correspond to the dash-dot line segment shown in Fig. [TJ Stable [unstable] fixed 
points and two-cycles are indicated by solid [dashed] curves. A supercritical period-doubling 
bifurcation occurs at fi = h^ijf) fa —0.2169; a border-collision bifurcation of the period- 
doubled solution occurs at /i = h^ 1 ^) ps —0.1937. When h^ 1 ^) < fx < there exists a 
chaotic attracting set. When fi > there is no local attractor. 
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admissible and unstable. Forward orbits are attracted to a chaotic solution created at the 
border-collision bifurcation at /i = 0. Unstable high period orbits are also generated at the 
border-collision bifurcation; these are of even period. To further illustrate the unfolding, 
Fig. [2] shows a bifurcation diagram corresponding to a one-parameter slice of Fig. [lJ-F. The 
chaotic attracting set collapses to the origin at [i = and to a two-cycle at fj, — h^ 1 ^). 

Two-cycles created in border-collision bifurcations at fi — consist of one point in each 
half-plane. In contrast, two-cycles created in period-doubling bifurcations along n = h\(n) 
lie entirely in the left half-plane. Along rj = h 2 (fj.) these two orbits coincide and are given 
explicitly by {0, /^(aO)}- The second iterate of (jHJ) can be transformed into a map 
of the form ([2]) along rj = h 2 (fi), see ( 1281) . The resulting left and right period- two maps 
have multipliers 1 — co/i 2 + o(/i 2 ) and —a^ + O(n) respectively. By Feigin's classification, 
four-cycles are born at r\ = /^(/u) exactly when > 1, i.e. in panels E and F of Fig. [TJ 
though these are not shown. 

4 Higher dimensions 

We now analyze the co dimension- two scenario in an arbitrary number of dimensions, i.e. x G 
R . Recall that for smooth systems the concept of dimension reduction by center manifold 
analysis allows bifurcations in systems of any number of dimensions to be transformed to their 
low dimensional normal forms. Unfortunately this technique cannot be applied here because 
([2]) is not differentiable in a neighborhood of the bifurcation. However we can compute a 
center manifold of the C k , left half-map. Assuming certain nondegeneracy conditions, we will 
show that the restriction of the left half-map to this center manifold yields a one-dimensional 
map that has the same form as the left half-map of (jSj) and satisfies the conditions fpl)-flm1) 
of Theorem [TJ Consequently the aspects of Theorem [TJ that incorporate only the left half- 
map extend immediately to the same scenario in higher dimensions. That is, there exists a 
curve f] = hi(fi) along which period-doubling bifurcations occur and a curve rj = h 2 (ji) along 
which the two-cycle created at the period-doubling bifurcation collides with the switching 
manifold. The drawback of this approach is that we fail to describe dynamical behavior of 
(J2J) that occurs on both sides of the switching manifold. 
Let us state our main result: 

Theorem 2. 

Consider the piecewise-C k , continuous map (TJP with (TJ) on R N and assume N > 2 and 
k > 4. Suppose that near (/i,^) = (0,0), Al^.ti) has an eigenvalue A(/x, 77) G R with an 
associated eigenvector, v(fi,rj). In addition, suppose 

i) A(0,0) = —1 is of algebraic multiplicity one and Al(0, 0) has no other eigenvalues on 
the unit circle, 

ii) f? T (0, 0)6(0, 0) 7^ (where g T is given by ^)), 
Hi) g(0,0) = l, 

iv) the first element of v(0,0) is nonzero, thus by scaling assume ejv(0,0) = 1, 
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v) det(/-A L (0, 0)^(0, 0))^0. 

Then, in the extended coordinate system (x; /J,,rj), there exists a C k ~ 1 three-dimensional cen- 
ter manifold, W c , for the left half-system that passes through the origin and is not tangent 
to the switching manifold at this point. Furthermore, f^\ Wc , described in the coordinate 

system (£; /L, 77) = (s; det 2 f_ b Al ^ L v)> given by the left half-map of with "hatted 
variables" and the conditions (fz))-/ fmj) of Theorem^ will be satisfied. 

//co / in Theorem^ (where Cq is computed for f^\ Wc ), then along a C k ~ 2 curve, 
f] = hi(fi), the unique fixed point of the left half-map near the origin undergoes a period- 
doubling bifurcation that is admissible when fi < 0, and along a C k ~ 2 curve, r) = h 2 (jx), the 
period- doubled solution collides with the switching manifold. These two curves intersect and 
are tangent at jx = fj = 0. Moreover, a unique two-cycle consisting of one point on each side 
of the switching manifold exists for small /i and rj. This cycle is admissible exactly when 

'det(I-A L )-JLdet(I+A L ) 



fx < and fj < sgn 



6st(I-A L A R ) 



(0,0) 



A proof is given in Appendix [Bl Conditions (TTT)— (TTTT|) are analogous to the first three 
conditions of Theorem [TJ Condition (jrvj) is a non-degeneracy condition that ensures the 
period-doubled cycles collide with the switching manifold in the generic manner. If this 
condition is not satisfied, a higher codimension situation that involves linear separation 
between and h 2 (fJ>) may arise. Finally condition (jvj) is necessary to ensure that non- 

degenerate, two-cycles are created at \i = 0. 

Theorem [2] essentially predicts the same basic bifurcation structure as Theorem HJ The 
relative position of hi and h 2 and the criticality of the period-doubling bifurcation are deter- 
mined by the sign of Cq. In one dimension, cq is a simple function of coefficients of nonlinear 
terms of the left half- map (1111) . To obtain the value of Cq in higher dimensions, one may 
derive an expression for the restriction of the left half-map to the center manifold, then apply 
the one dimensional result. This is done for an example below. An explicit expression for cq 
in N dimensions is probably too complicated to be of practical use. 

To illustrate Theorem [2] consider the following map on M. 2 using x = (s,y): 

1 1 

s = -9^- ~ s + y , 

, 1 i, , 3 1, < ie > 



This map is easily rewritten in the piecewise form (j2J) with ([3]) and as we will show ( 1161) 
satisfies the assumptions of Theorem [2j When \i = the origin is a fixed point of (fl6|) which 
lies on the switching manifold, s = 0. Consequently the r^-axis corresponds to a locus of 
border-collision bifurcations. 

Fig. [3] shows a two-dimensional bifurcation diagram of (1161) near fj, = 77 = 0. Period- 
doubling bifurcations occur along rj = h\{ix) and the period-doubled solutions undergo 
border-collision along 77 = h 2 (ti)] these curves are given explicitly by 

\ 48 (17) 
hz(fi) = --fi , 
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(actually (rj) = —6r/ — fr/ 2 ). The curves are tangent at \i = r\ = as predicted by 
Theorem [2J Near this point the dynamics resemble Fig. [Tp although there is an additional 
unstable two-cycle present that does not undergo bifurcation at /j = r\ — 0. This two-cycle 
collides with the stable two-cycle born at fi — for small rj < along a locus of saddle- 
node bifurcations (shown red in Fig. [3]). The saddle- node locus emanates from the r?-axis 
at rj = ~ 4 +^° rj _o.093 which is where the stable two-cycle loses stability. Also a cusp 
bifurcation of period-two orbits occurs at (/-t, 77) = (—j^, — §)■ 

A second simultaneous border-collision and period-doubling bifurcation occurs on the 
77- axis at 77 = — |. As predicted by Theorem [2] a locus of period-doubling bifurcations and 
border-collision bifurcations of the period-doubled solution emanate from the codimension- 
two point and are tangent to one another here. Locally the bifurcation set resembles Fig. [T]A 



'1 




-0.5 0.5 1 1.5 

Figure 3: A bifurcation set of the two-dimensional, piecewise-C°° map (fl6l) . Border-collision 
bifurcations of fixed points and two-cycles are indicated by black curves which are dashed 
when no bifurcation occurs in a topological sense. Blue curves correspond to period-doubling 
bifurcations of a fixed point; red curves correspond to saddle- node bifurcations of a two-cycle. 
As in Fig. [H schematics showing fixed points and two-cycles are included. Dynamics local 
to the codimension-two points (/j, 77) = (0,0) and (0, — |) are described by Theorem [2j 



We now perform calculations to verify the details of Theorem [2] for the map (TT6T) near 
fj, — Tj — 0. Here £ T (0, 0) = [1, 1], 6(0, 0) = [— |, 0] T and det(I - A L (0, 0)) = 1, thus fx = -/j 
and therefore hi(fi) corresponds to admissible solutions when /i > 0, matching Fig. [31 Also, 

det(J + A L ) =| and det(J - A L (0, 0)A R (0, 0)) = |, hence by the last statement of 

Theorem [2J two-cycles born on /i 2 with points on each side of the switching manifold are 
admissible below /i 2 , in agreement with Fig. [3j 
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Using a series expansion (see also the proof of Theorem [2]), one finds that the left half- map 
of f)16p restricted to the center manifold is given by 



fi - s + 



1 „ 2 

-8 



7]S+ -fi 



2flT] 



:s 3 + 



Here the coefficient —(jpr- + p) , given in Theorem[[J is equal to ~. Also c = — ^. Thus 



(0,0) 



hi(fi) = — + 0(/i 2 ) and h 2 (fi) — hi(fi) = -hfi 2 + 0(fi 3 ) matching the explicit results given 
above. 



5 Conclusions 

A two-parameter bifurcation diagram provides a useful and concise picture of the behavior 
of a dynamical system near a codimension-two point. Loci or curves in such diagrams corre- 
spond to codimension-one phenomena. The intersection of two or more bifurcation curves at 
a codimension-two point occurs in a manner that is generally determined by the correspond- 
ing codimension-one bifurcations. We have shown that for a general piecewise-smooth, con- 
tinuous map, loci of border-collision and period-doubling bifurcations of a single fixed point 
generically intersect non-tangentially (unlike the analogous intersection of a border-collision 
bifurcation and a saddle-node bifurcation [IB])- Moreover the period-doubled solution un- 
dergoes border-collision along a locus that emanates from the codimension-two point and is 
tangent to the period-doubling curve here. 

We have determined all local dynamical phenomena that occurs when this scenario arises 
generically in a one-dimensional system. In particular we have shown exactly when chaotic 
dynamics may be generated. In higher dimensions border-collision bifurcations near the 
codimension-two point may in addition generate quasiperiodic solutions such as invariant 
circles or other more complicated orbits. 



A Proof of Theorem H 

The theorem is proved in five steps. First, the function h\ is calculated by finding where 
x*^ has an associated multiplier of —1. Second, h 2 is calculated from an expression for 
f( L ' . Third, an explicit computation of the piecewise-smooth, second iterate map allows 
the border-collision bifurcation of the two-cycle that occurs along 77 = h 2 (fi) to be classified. 
Finally, in steps four and five we prove parts (jlj) and (J5J) respectively, 
since © is C k and k > 3, we can write 

a L (fi, 77) = -1 + aifi + 77 + a 3 fi 2 + ctAfir} + a 5 r/ 2 + o(2) , 
b(fL, V ) = l + /3 lf jL + f3 2 ri + 0(2), 

P((*,V) = 7o + 7iA* + I2V + o(l) , 
q(fi,r)) = 5 + o(0) , 

where we have used hypotheses of the theorem to simplify the coefficients. 
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Step 1: Compute the function hi which corresponds to the existence of a fixed point of 
the left half-map with multiplier — 1. 

By the implicit function theorem the left half-map has the fixed point 

(n, rj) = -fi + hfi 2 + k 2 /iri + k 3 r] 2 + 0(3) , (19) 

which is a C k function of \x and rj and locally satisfies f^ L \x*^ L \fi, rj); /i, rj) — x*( L \[x, rj) = 0. 
By a second order expansion, it is determined that the scalar coefficients, fcj, have the values 

, A . oli ,7o 

fci = y + T + T' 

= * + i < 20 > 

2 4' 
A: 3 = 0. 

The multiplier associated with x*^ L \fj,,rj) is found by substituting ( fT9l into D x f( L > (x; fi, rj) 
(see (ED). We obtain 

D x f( L \x< L \v,v);tJ>,v) = -1 + ("i + 7o)/i + »7 + (as + 2A;i7o + 7i + ^f) M 2 

+ (o 4 + 2A; 27o + 7 2 )/xr/ + c^ 2 + o(2) . (21) 

The implicit function theorem implies that 

hM = -(^ + 7o )// + /i/i 2 + o(// 2 ) , (22) 

is a C^ -1 function that locally satisfies D a ./^ L ^(x*^ L ^(/i, hi(fi)); /i, hx(fj,)) + 1 = 0, where /i is 
found by substituting (J22D into (121)1 . 

/i = - + 2£;i7o + 7i + + (a 4 + 2fc 2 7o + 72) («i + 7o) - ot 5 (ati + 7o) 2 • (23) 

When Co 7^ 0, the existence of period-doubling bifurcations along r\ = hi(fi) for small fi is 
verified by checking the three conditions of the standard theorem, see for instance [33J: 



i) (singularity) by construction, D x p L \x*^ L '{ii 1 hi(fi)); fi, hi(fi)) 

ii) (transversahty) + 2^) | _ = 2 + 0(/i) + 



iii) (non- degeneracy) I ^ y 



1 ( d 2 fW \ , 1 a 3 /< L > 

9x 2 / 3 9x 3 



= 2c + O(fi)^0. 



ri=hi(n) 



Consequently we have proved ([T]) of the theorem. 

Step 2: Compute the function h 2 which corresponds to the existence of a period-two 
orbit of the left half-map at x = 0. 
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To determine h 2 , we compute the C k function /^ 2 (0; fi, rj) = /^ L ^(/i6(/i, 77); /i, rj) which may 
be written as f^ L ' (0; /i, rj) = fig(fi, rj) where 

g(fi, 1]) = 6(//, rj) (l + a L (/i, rj) + rj)p{ji, rj) + /i 2 & 2 (/x, r))q(/i, 77)) + o(2) , (24) 

is C fc_1 . Substitution of ([181) into (Ell) produces 

flr(/i, 77) = («! + 7o)/j + 77 + (aiA + a 3 + 2/3i7 + 7i + 5 oj> 2 

+ (A + a x (5 2 + a 4 + 2/?2 7o + 72)/^ + (A + «s)?7 2 + o(2) . (25) 

The implicit function theorem implies that 

^2 (a*) = + 7o)a* + kfJ 2 + o(/j 2 ) , (26) 

is a C fc_1 function that locally satisfies g(ji, h^ji)) = 0. By substituting (1261) into (1231) we 
deduce 

Z 2 = -(«i/5i + a 3 + 2/?i7o + 71 + 5 ) + (A + ai/5 2 + a 4 + 2/5 2 7o + 72)(«i + 7o) 

- (& + a 5 )(ai+7o) 2 • (27) 

Subtracting f|23]) from (127]) produces (after algebraic simplification) 

, 7 _ To + ^0 
«2 - k - — , 

which proves (1I4"1) in the statement of the theorem. 

Step 3: Determine all two-cycles to verify the phase portraits in Fig. [3 
For small, fixed ji < 0, the two-cycle generated in a period-doubling bifurcation at rj = hi(fi) 
undergoes a border-collision bifurcation when it collides with the switching manifold at 
77 = /i 2 (/i). The stability and relative admissibility of two-cycles emanating from this border- 
collision bifurcation is found by determining a map of the form (1131) that describes the 
bifurcation. Such a map may be obtained by computing the second iterate of (jHJ) and 
replacing 77 with the new parameter 

77 = 77 - h 2 (n) , 

which controls the border-collision. Near 77 = 0, when \i < period- two orbits are guaranteed 
to have one negative point, so we compute o and o which leads to 



2 J (/i + 0(/i 2 ))77+(l-coA4 2 + o(/j 2 ))x + 0(|x,77| 2 ), x<0 

1 i ^ + o(n 2 ))fj + (-a^ + 0(ji))x + 0(\ X ,fj\ 2 ), x>0 



A fixed point, x*^ LL \ of the left half-map of (1281) corresponds to a two-cycle of (jHJ) with 
both points negative- valued, when admissible. By (|28l) . if Co < [cq > 0], then 

x *(ll) ig 

unstable and admissible when 77 < [stable and admissible when 77 > 0]. Similarly, a fixed 
point, x*( RL ^ of the right half-map of (1281) corresponds to a two-cycle of (jHJ) comprised of two 
points with different signs, when admissible. If \o,q R \ < 1, then x*^ RL ^ is stable, otherwise 
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it is unstable. Furthermore if Oq < — 1 [a^ > — 1], then X *( RL ) is admissible when fj > 
[77 < 0]. Lastly, since ag 7^ —1, the fixed point, x*^ L \ is unique for small // and 77. 
Therefore the two-cycle created at \x = collides with the switching manifold at 77 = h 2 {fi). 
These statements verify all two-cycles shown in Fig. [TJ 

Step 4: STicw £/iai /, (TJj] ; has no n-cycles with n >3, verifying part Oj. 
Except when a = our proof is founded on the knowledge that a one- dimensional map 
cannot have an n-cycle with n > 3 contained in an interval on which the map is monotone. 
For the special = we used several additional logical arguments including the 

contraction mapping theorem. 

Case I: a ( R) < 0. 

In this case there exists 5 > such that V/i, 77 E [—5,5], f is decreasing on [—5,5]. Conse- 
quently / has no n-cycles with n > 3 on [—5, 5]. 
Case II: a { R) = 0. 

First we construct an interval containing the origin on which / is forward invariant. Since 
the two components of / are differentiable, there exists 5 > such that V/x, r],x G [—5,5], 



d f(R) 



dx 



Qf(L) 

+ 1 



dx 



\b(fM, V )-l\< 1 



3 



Let 5 = |5, let I = [—45,85] and assume 6 [—5,5]. Then /(/) C / because for any 
x £ I, we may assume x > when calculating a lower bound for f(x): f(x) > fib — 
|x > — 15 — |5 = —45, and we may assume x < when calculating an upper bound: 
f(x) < fib - fx < |5 + |45 < 85. 

If (i > then the image of any x < under / is positive and since the right half-map is 
contracting and the slope of the left half-map is near —1, f 2 is a contraction on J, indeed 
\D x f 2 (x)\ < 1 3 < 1. By the contraction mapping theorem, iterates of / 2 approach a unique 
fixed point, which must be x*^ R \ therefore / cannot have an n-cycle with n > 3 on J\f (|12p . 

So we now assume fi < 0. Recall that when 77 = /i2(/i), / ( - L - ) (0) = 0, i.e. /(0) = (0), 
and when 77 > h 2 (fi), /(0) > (0). Therefore if 77 > h 2 ((i), see Fig. 01-A, the interval 

(0), 0] is forward invariant. Since / is monotone on the interval 

J =[/(()), 0], (29) 

it contains no n-cycles with n > 3 here. On the other hand if 77 < ^(a*) as in Fig. H]-B, 
the forward orbit of any point in [f^ 1 (0),0] is either contained in this interval (where / 
is monotone) or enters [/(0),/( L ) (0)] and approaches the stable two-cycle (this may be 
shown formally by the contraction mapping theorem). In any case, an n-cycle of / with 
n > 3 cannot include a point in J = [/(0), 0]. 

Suppose for a contradiction / has an n-cycle with n > 3, call it {x . . . x n _i} C [—5, 5] \ J. 
We now show that the restriction of f 2 to this n-cycle is a contraction giving a contradiction 
according to the contraction mapping theorem. Specifically we will argue that for any i and 
j, 

\f( Xj )-f( Xi )\<-\ Xj - Xi \. (30) 
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There are several cases to consider. For example, suppose that Xj > 0, Xi < and f 2 



< 



f\x t ). Then f( Xj ) > /(0) - \x 3 thus f 2 {x 3 ) > f(0) + ±/(0) 

1 

7'^ A 



thus /2( Xi )< / ( ) + I/(0) 



9 x r 



Also < /(0) - f a* 



Subtracting these inequalities leads to (13"Uj) . The remaining 



cases can be shown similarly. 

u > h 2 {\k) 




Xn ■ l 



\ 




/' L ' J (0)X\ 






/(»K ^ 






J 


K 





Figure 4: Schematic diagrams illustrating the map /, ([8]), when /i < and 77 > /t-2 (/^) in 
panel A and r\ < h 2 (n) in panel B. When 77 = /i 2 (/i), 1 (0) = /(0). 



Case III: of' > 0. 

Case Ilia: First suppose that /i > 0. There exists Si > such that whenever /z G [0, <5i] and 
77 G [— £1, / is decreasing on [—Si, 0] and increasing on [0, Si}. Since /i > 0, Vx G [— 
f(x) > 0. Thus any periodic solution of / in [— Si, Si] must lie entirely in [0, Si]. But / is 
increasing on this interval thus has no n-cycles with n > 3 here. 

Case Illb: Now suppose \i < and 77 > h 2 (/i). There exists S 2 > such that whenever 
/i G [—S 2 ,0) and 77 G [/12 (aO ? ^2] 5 / is decreasing on [— 5 2 ,0] and increasing on [0, #2], see 
Fig. HJ-A. Suppose for a contradiction, / has an n-cycle with n > 3 on [—S 2 ,S 2 ]. Such 
an n-cycle must enter both [— <5 2 ,0) and (0,o" 2 ], and so it includes a point x G (0, S 2 ] with 
f(x) G [S 2 ,0). But f(x) lies in J, ([21]), and since 77 > /i 2 (m)> / 2 (0) < 0, hence /(J) C J. 
Thus the forward orbit of x cannot return to (0, S 2 ] which contradicts the assumption that 
x belongs to an n-cycle. 

Case IIIc: Finally suppose /1 < 0, 77 < h 2 ([i) and < a Q < 1. Then there exists £3 > 
such that whenever jj, G [—S3, 0), 77 G [—S3, h 2 (/i)) and x G [—£3, S3], 

(ii) , (i?) 

^< fW(x) <±±^°-, 

2 v ; - 2 

where the particular non-symmetric -dependent bounds on the slopes have been chosen 
to simplify the subsequent analysis. Suppose for a contradiction / has an n-cycle with 
n > 3 in [—S3, S3]. As before, such an 77,-cycle must include a point x G (0, £3] with f(x) G 
[— S 3 ,0). Again f(x) G J, but here 77 < h 2 (fi), thus the forward orbit of x may return 
to the right of the origin but must enter the interval K = [0,/ 2 (0)], see Fig. 0]-B. Notice 
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/ 3 (0) < /(O) + ^^f 2 (0). Then by evaluating the first-order expansion of about /(0) 

(1 («) \ 1 i (R) 
1 + I 1 ^) 1± | L / 2 (0) > 0, therefore / 2 (K) C K. 
l+a y 

/ / , (H)\ u (H) 

Then Vy G if, |/ 2 (y)\ < 1 + l 1 ^, i±fB_ < 1. Thus / 2 : K -> if is a contraction 

\ 1+a / 

mapping and therefore iterates of / that enter if cannot belong to an n-cycle with n > 3. 
Finally let 5 = min(5 1 , 5 2 , S 3 ), then the result, (|S]), is proved. 




Figure 5: The second iterate of (j2j), / 2 , near a; = when p < 0, 77 < ^(a 4 ) an d aj^ > 1. The 

f 

slope of this map is approximately unity left of the origin and approximately — % right of 
the origin. As a consequence, see the text, the map exhibits chaos in T = (f(2p)), 2p 



Step 5: Prove the existence of chaos when /x < 0, rj < /^(aO an d % > 1. 
We first construct a trapping set, T, for / 2 . Let p(p,n) = / 2 (0; p, 77). Since p < 0, by (T28j) . 
p(p, 77) = 0(2). Since 77 < h 2 (p), p > 0. Therefore 

/(2p) = /(0) + (2a^ } + 0(l))p + 0(p 2 ) < . 

Consequently, 

/ 2 (2p) = (l-2a^+0(l))p + 0(p 2 ), (31) 
f^-\f(2p)) = (-2af + 0(l))p + 0(p 2 ). (32) 



Let 



T 



> w l (/(2p)),2p 



It follows, see Fig. El that Vx G T, / 2 (2p) < f 2 {x) 



< 



But p,/ 2 (2p) G int(T), thus 



/ 2 (T) C int(T), i.e. T is a trapping set for / 2 . Furthermore, fli^o / 2? (-0> i s an attracting 
set for f 2 . 
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Let M G Z with M > 2a ( R] + 2. We now show that Va; G T, 3j G Z with < j < M 
such that f 2 ^(x) > 0. Suppose otherwise. Then whenever < j < M, / 2j (x) < and thus 

f 2j (x) = (j + 0(l))p + (1 + 0(l))x + 0(x 2 ) . (33) 

Since f 2 is increasing on T_, it suffices to consider x = 1 (/(2p)). By combining ( 1321) 
and fl33l) we find 



J 2 ^- 1 )^) > (24 /?) + l + 0(l))p+(l + 0(l))(-2a^ ) +0(l))p + 0(p 2 ) 
= (1 + O(l))p + O(p 2 )>0, 

which is a contradiction. 

Now consider the map f 2M : T — > T. Note that J 2 has one critical point on T, namely 
0. Thus y G T is a critical point of f 2M if f 2 i{y) = for some < j < M. Consequently 
f 2M has at most 2 M — 1 critical points, yi, and f 2M is different iable on T except at each yj. 
For all x G T, x ^ yi, 

M-l 

f M \x) = n f\f 23 ^)) ■ 

3=0 

Let m > 1 be the number of iterates, f 2 ^(x), that are positive. Since left of zero the slope 
of f 2 is near 1 and right of zero the slope is near —a^ , we have 

/ 2M >) 



^ )m + 0(l)>i^>l 



for any x ^ yi, for sufficiently small /i, r/. Consequently J 2A/ : T — > T is piecewise expanding 
[3l]. By the theorem of Li and Yorke [35], f 2M is chaotic on T. Thus / is chaotic on 
TU/(T), □ 



B Proof of Theorem 2 



The theorem is proved in three steps. In the first two steps a center manifold, W c , for the 
left half-map, f^ L \ is constructed and the lowest order terms of the restriction of to W c 
are calculated. In the third step the two-cycle that has points in each half-plane is computed 
explicitly. 

Step 1: Compute W c . 



It is convenient to extend the map to 
phase space, 



X 



by including the parameters in the extended 





x' 




' f L \x;fi,rj) 


F = 


\£ 








rf 




V 



Then 



DF(0;0,0) 



^4l(0,0) 


6(0,0) 








1 











1 
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has a three-dimensional center space 



where 
Notice 





' v{0,0) ' 








E c = span | 







1 
















tp=(I 


-A L (0,0))- 


x 6(0,0). 



£ T (0, 0)6(0,0) 



7^0. 



(34) 
(35) 



det(J-A L (0,0)) 

Since ejv(0, 0) = 1 ^ 0, the center manifold, W c , may be expressed locally in terms of s, ji 
and rj. By the center manifold theorem, there exists a C fc_1 function, if, such that on W c 

x = H(s; /i, rj) = sv(0, 0) + fiC + 0{2) (36) 

where 

{=ip-(ej<p)v(0,0). (37) 

Step 2: Determine an expression for on W c and verify conditions /fzj)-/ fmj) of The- 
oremUl 
On W c 



s' = e[f^{H{s-^r,)-^r,) 

= ejb(n, r))n + e[A L (n, r))H(s\ ji, rj) + e]g {L) (H(s; n, 77); //, 77) 



(3? 



where denotes all terms of f^ that are nonlinear in x. With the given "hatted" variables 
satisfies condition (TTT1) of Theorem [T] because 



ds' 

= eJb(0,0) + elA L (0,0)C 

Oji (0;0,0) 

= ejb(0, 0) + ejA L (0, V - (eJ<p)ejA L {0, 0)v(0, 0) (by (J37D) 

= e^ + (e^)e^(0,0) (by dMD) 

= 2e]"y?, 

which with ji = 2eJ \pfi, in view of (j3"5|) . produces 

5/2 (0;0,0) 

A similar verification of (J]} and flm|) of Theorem [1] is now given. By scaling v, we may assume 
that for small rj, ejv(0, rj) = 1. Since fl36l) lacks a linear rj term, 



and consequently 



H(s ] 0,r i )=v(0,r i )s + O(s 2 ) , 

U = o = ejA L {Q,rj)v{Q, rj)s + 0{s 2 ) 
= X(0,rj)s + O(s 2 ) . 
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Thus 



ds 



(0;0,0) 



A(0,0) 



d 2 s' 



drjds (o ; o,o) dr] 



(0,0) 



-1 
1 . 



as required. 

Step 3: Compute the two-cycle that has points in each half-plane to verify the final 
statement of the theorem. 
We have 



fM(fW(x; fx, 77); fi, ti) = {I + A L (0, V ))b(0, n)n + A L (0, r))A R (0, rf)x + 0(\(i, x\ 



(39) 



By applying the implicit function theorem to fl39|) in view of condition (jyj) of the theorem 
we see that the two-cycle with points in each half-plane exists and is unique for small \x and 
r). When admissible, the point of this cycle in the right half-plane is given by a C k function 



x 



*(RL), 



Ou, 77) = (I — A L A R )-\I + A L )b 11 + O(^) . 



Then 



Q S *(RL) 



At=0 



ej(l - A L A R )-\I + A L )b 
e]adj(/- A L A R )(I + A L )b 



fj,=0 



det(7 — A L A R ) 



Since I — AlAr = (I + Al)(I — A R ) — Al + A R and and A# are identical in their last 
N — 1 columns, adj(J — AlA r ) and adj((J + Al)(7 — Ar)) share the same first row. Thus 



Q S *(RL) 



e[adj(/ - A fl )adj(7 + A L )(/ + A L )6 



/x=0 



det(J — AlAr) 
det(J + Al) 



-g T b 



det(J — AlAr 

jdet(J + A L ) 
det(J — AlAr) ^ 



fi=0 



V + 0(rf 



(0,0) 



By Theorem [T], there exists a C k 2 curve, 97 = /2.2(A)) along which, s*( RL ) = 0. Furthermore, 
when fi < 0, the two-cycle is admissible along this curve. If 7] = 77 — /12(A)) then 



(A^) 



det(/-A L )|^det(/ + A L 



2det(/-A L A R ) 
which confirms the final statement of the theorem. 



Ar] + 0(3) , 



(0,0) 



□ 
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